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Abstract 

I calculate the top quark forward-backward asymmetry at the Tevatron in both the lab¬ 
oratory frame and the tt rest frame. I show that soft-gluon corrections are the dominant 
contribution to the asymmetry and closely approximate exact results through next-to- 
next-to-leading order (NNLO). I present a calculation of the asymmetry including approxi¬ 
mate next-to-next-to-next-to-leading-order (N 3 LO) soft-gluon contributions from next-to- 
next-to-leading-logarithm (NNLL) resummation as well as electroweak corrections. Thus 
approximate N 3 LO (aN 3 LO) results are obtained, which significantly enhance and improve 
previous NNLO results. The theoretical aN 3 LO result for the top quark forward-backward 
asymmetry at the Tevatron in the laboratory frame is (6.8±0.3)%, and in the tt rest frame 
it is (10.0 ± 0.6)% which is in excellent agreement with recent Tevatron data. 


1 Introduction 

Top quark production has been studied extensively at the Tevatron and the LHC. The soft-gluon 
corrections are now known approximately to next-to-next-to-next-to-leading-order (N 3 LO) [lj 
[2], based on the next-to-next-to-leacling-logarithm (NNLL) two-loop soft-gluon resummation 
formalism of [3], which allows the calculation of approximate N 3 LO (aN 3 LO) results for the total 
cross section [lj and the transverse momentum and rapidity distributions [2]. The agreement 
between theory and experimental results from the LHC and the Tevatron is excellent for all 
available energies for both the total and differential cross sections. 

A quantity of particular interest is the top quark forward-backward asymmetry, Afb ; which 
is closely related to the rapidity distribution of the top quark. The definition of the asymmetry 
is 

i FB = An > °> - Ai < ^, (!.!) 

a(y t > 0 ) + a(y t < 0 ) a 

where a denotes cross sections and yt is the top-quark rapidity in a given Lorentz frame, which 
may be the pp laboratory frame or the tt rest frame (in the latter frame the asymmetry is 
equivalently defined in terms of Ay = y t — yi, with yi the antitop rapidity; see Ref. [5] f° r more 
details on various definitions and further references). At a pp collider such as the Tevatron, 
charge conjugation invariance in QCD implies that a{yi > 0) = a(y t < 0), so the forward- 
backward asymmetry is equivalent to the charge asymmetry. The tt rest frame asymmetry is 
larger than the laboratory frame asymmetry by about fifty percent due to kinematics. 

Early measurements of Afb were performed by CDF [5] and DO [6] and returned values 
substantially larger than standard model predictions. Recent measurements of Afb hi the tt 
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rest frame with more data and precision have been reported in [0, 8]. Thus, it is important 
to perform the most accurate calculation for the asymmetry in order to have confidence in the 
theoretical prediction while looking for hints of new physics. 

The QCD corrections to the asymmetry were calculated with next-to-leading-logarithm 
(NLL) resummation of soft-gluon contributions in the moment-space formalism in [9]. This 
was later extended to NNLL accuracy in [TO] where approximate next-to-next-to-leading-order 
(NNLO) results were obtained. Results at NNLL were also derived using Soft-Collinear Effective 
Theory (SCET) in [IT]. An approach using the principle of maximum conformality was used 
in ra to calculate the asymmetry. Recently the complete NNLO QCD contribution to the 
asymmetry has appeared in m- 

A variety of resummation formalisms and methods have been proposed in the literature 
for top quark production over the past twenty years (see the discussions in (3[ EH [15, 16] and 
references therein). These differences include whether resummation is performed for the total 
cross section or for double-differential cross sections (i.e. absolute vs partonic threshold); the 
inclusion and motivation for various subleading terms and whether or not they are warranted; 
the use of damping factors away from threshold; whether resummation is performed in moment- 
space perturbative QCD vs SCET resummation; the use of prescriptions (and their validity) 
in resummed results vs performing prescription-independent fixed-order expansions at NNLO 
or N 3 LO; the choice of kinematics, single-particle-inclusive (1PI) vs pair-invariant-mass (PIM), 
for double-differential resummations, and various choices for writing the analytical terms via 
partonic threshold relations of kinematical variables; etc. More details and references may 
be found in J3[ [EH ESI EE]- In most recent formalisms (see Refs. [3] EH ESI ESI EH EH]) the 
soft-gluon corrections were found to be large and dominant. 

Most resummation formalisms for top production use methods applicable only to calcula¬ 
tions of total cross sections. Only two formalisms have been developed at NNLL for double- 
differential cross sections and applied to the asymmetry, one using moment-space QCD resum¬ 
mation Rij and the other using SCET m- The differences between the moment-space QCD 
and the SCET resummation approaches in top production were discussed in the review paper 
[15], where results for the asymmetry in different formalisms were shown. In the formalism of 
m the approximate NNLO contributions to the asymmetry were found to be significant while 
in the approach of mi they were rather small. This discrepancy was not unexpected since the 
overall contributions from soft-gluon resummation to the total cross section and differential dis¬ 
tributions were also quite different in the two approaches. There are several differences between 
the two resummations, as described in [15], but the numerically prominent one is the inclusion 
of certain terms in m that arise in the SCET approximation and are not strictly of soft gluon 
origin. Those SCET terms severely depress the results in m relative to HD]. As we will show 
in Section 2, the result with large approximate NNLO corrections in [10] precisely predicted 
the exact NNLO asymmetry; the results in [TT] are much smaller, with tiny corrections beyond 
next-to-leading order (NLO). 

The fact that the resummation formalism used in [3j EE] predicts the exact NNLO cross 
section better than all other approaches, to better than one per mil accuracy for the central 
value at the Tevatron, as shown in [1(5], had already strongly indicated that the results for 
the asymmetry in [TO] would be very close to the exact NNLO result before the latter was 
known. This is indeed the case as we discuss in the next section. We note that there are strong 
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theoretical reasons for the choices we have made in onr approach. Formally, the resummation 
uses renormalization group methods to fully determine terms at a given logarithmic accuracy 
(NNLL in this case). Numerically, subleading terms have to be treated carefully to approximate 
and predict the corrections correctly. At NLO the total cross sections as well as the transverse 
momentum and rapidity distributions at all Tevatron and LHC energies were reproduced nearly 
exactly with our method, ft was also shown that with our method the differences between 
approximate NNLO results in 1PI and PIM kinematics are minimized, which indicated that 
any missing terms were numerically small [T9] , This fact together with the correct description 
at NLO was a very strong indicator that our NNLO predictions would be successful and, not 
surprisingly, this has indeed been the case, to very high accuracy, which validates our method 
and the theoretical arguments in support of it. 

Parenthetically we note that it has been known for a while that the NLO expansion of the 
results from soft-gluon resummation for the asymmetry reproduces the complete NLO result 
rather well (this was noted already in [9] where NLL resummation was used). The extension 
to NNLL in US] has a similar effect in the NNLO expansion, in predicting the exact NNLO 
results of [13] very well. 

The authors of [13] also found a significant NNLO contribution, as expected from [10] . 
However, the authors of [13] stated erroneously that the effects of soft-gluon contributions are 
negligible. As discussed above, and as will be shown in more detail in the next section, this is 
factually incorrect. The soft-gluon contributions have been known to be large and dominant in 
the total and differential top production cross sections for at least twenty years (for references 
see the review in [15]), and in the asymmetry in particular at least since 2011 [10]. In fact the 
NNLO enhancements to the asymmetry found in [13] are very close to the predictions in [ID] . 
The accuracy, precision, reliability, and robustness of the resummation formalism that we use 
have been amply demonstrated and discussed in [3, IB] for total cross sections and differential 
distributions. 

We also note that while our resummation formalism uses the soft-gluon eikonal approxima¬ 
tion for partonic threshold, the integrations are over all allowable gluon energies. Furthermore, 
partonic threshold is more general than absolute production threshold, i.e. the top quark is 
not necessarily produced at rest but may have arbitrary transverse momentum and rapidity as 
allowed by the kinematics. As a consequence, the effects of these approximate contributions to 
the top-quark transverse momentum (px) distribution increase with px at large px since that 
is a region closer to partonic threshold (see e.g. Fig. 3 in the first paper of Ref. [16]). Thus, 
we find the arguments made in [13J regarding the contributions of soft gluons to the transverse 
momentum distribution of the top-antitop pair to be invalid or irrelevant in our approach. 

In this paper we will improve on the NNLO results for the asymmetry by further calcu¬ 
lating N 3 LO soft-gluon corrections, and including electroweak corrections [H [20] [2l] [22], thus 
producing an aN 3 LO estimate. We present the results for the aN 3 LO asymmetry in both the 
laboratory frame and the tt rest frame. 
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2 The top-quark forward-backward asymmetry at the 
Tevatron 


We write the perturbative QCD expansion for the top-antitop pair cross section as 


a = 5> s 2 +V n > 

n=0 


( 2 . 1 ) 


with a s the strong coupling. At leading order the cross section is of order ex 2 , the NLO cor¬ 
rections ed 1 ) are of order a 3 , etc. We include the complete nd) an d a ( 2 ) corrections and the 
approximate cd 3 -* corrections from NNLL resummation. 

At each order in cx s , the cross section includes plus distributions of the form [ln A: (s 4 /m 2 )/s 4 ] + 
where S 4 = s+t+u— 2 m 2 , with s, t, u the usual kinematical variables and m t the top-quark mass. 
The variable S 4 measures distance from partonic threshold and, for the nth-order corrections, 
the power of the logarithm, k, in the plus distributions can range from the leading value of 
2 n — 1 down to the lowest value of 0. Thus, at NLO the leading value of k is 1, at NNLO it is 
3, and at N 3 LO it is 5. With NNLL resummation one can derive the coefficients of all powers 
of the logarithms at NLO and NNLO, but not at N 3 LO. We use the inverse transform of the 
moment-space logarithms of the Mcllin moment to the specified order, including constant Q 
terms (see also discussion in P31). As discussed in p], at N 3 LO the coefficients of the powers 
of the logarithms with k = 5,4, 3,2 can be fully determined, but those with k — 1,0 are only 
partially known and include factorization and renormalization scale terms as well as terms 
arising from the inversion from moment space back to momentum space (see Ref. |23j for 
explicit expressions). The inversion terms, involving Q constants, are known to dominate the 
subleading powers of the logarithms, and the impact of the remaining subleading terms is much 
smaller than the scale variation [ 1 ]. 

The top quark forward-backward asymmetry is then written as 


Acr EW + cx 3 And) + ad Aa (2) + a 5 s An (3 ) + • • • 
a 2 n(°) + a 3 nd) + a f nd) + a b s nd) + • • • 


where again Acd n ) denotes the QCD asymmetric terms as in Eq. (11.11) . i.e. the difference of 
nth-order QCD corrections integrated over positive and negative rapidities, and Acr EW denotes 
the asymmetric term from electroweak corrections. 

One may further reexpand the above ratio in the definition of Afb in powers of a s , thus 
obtaining the expression 


^4fb 


A(7 EW 

-b OL 

a 2 cd°) 

Aa EW r 

+ (^M) 3 ■ 

-f- 


And) 
s <d°) 


(n (1) ) 2 


Acr EW crd) 

-b OL 

cx s (n (°)) 2 s 


- cd°) ad) 


+ 


And) 
n(°) 
And) 


n 


( 0 ) 


And) a (D- 
(n (°)) 2 

And) An«(nd )) 2 

(n (°)) 2 + (n (°)) 3 


And) a (2)- 
(n(°)) 2 

(2.3) 


Explicit formulas and numerical values for the quantities And) anc } a w ma y be found in 
Refs. jH 0, :T3] [23]. In the following we will provide results for both definitions, Eq. (I2.2j) and 
Eq. (12.31) . in both the laboratory and tt rest frames. 
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In the Standard Model the forward-backward asymmetry is due mainly to QCD effects. 
At leading order in perturbative QCD, the production channels qq —> tt and gg —> tt are 
symmetric in rapidity, thus App vanishes at that order. Furthermore, the gg channel remains 
symmetric at all orders in perturbative QCD. However an asymmetry arises in the qq channel 
(as well as a much smaller contribution from the the qg channel) starting at NLO. Therefore by 
applying soft-gluon resummation we expect the gg channel to remain symmetric, but we have 
contributions to the asymmetry from higher orders in the qq channel. 

2.1 pp laboratory frame 

We first provide results for the asymmetry in the laboratory frame. We use MSTW2008 [24] 
parton distribution functions and a top quark mass mt = 173.3 GeV in the calculations. The 
central values refer to the choice of scale /i = m t , and theoretical errors are estimated with 
scale variation by a factor of two. 

Using the NNLO approximate rapidity distributions it was found in ra that the approxi¬ 
mate NNLO top quark forward-backward asymmetry in the laboratory frame at the Tevatron 
(from QCD effects only) is 5.2%, using Eq. (12.21) . This is a large increase of thirty percent over 
the NLO asymmetry (and as we will see later, we find essentially the same percent increase in 
the tt rest frame). We also note that this NNLO / NLO ratio is very similar to the increase 
that was found in [13], as we also discuss in the next subsection. We confirm the results of 
m in this paper and we add the aN 3 LO QCD corrections as well as electroweak corrections. 
Furthermore, we calculate the corresponding numbers using Eq. (12.3[) . 

At aN 3 LO we find that the asymmetry from QCD alone is (5.61 q; 4)% using Eq. (12. 2\ . and 
(6.0 ± 0.1)% using Eq. (12.31) . The uncertainties indicated in these numbers are estimated by 
varying the scale p between m t /2 and 2 m t . If we include the electroweak corrections as well as 
the aN 3 LO QCD corrections we find for the asymmetry the results (6.4 ±qq)% using Eq. (12.2|) . 
and (6.8 ±0.3)% using Eq. (12.3p . The aN 3 LO results with and without electroweak corrections 
are tabulated in Table 1. 

We also note that the differences between using Eq. (12.21) and Eq. (12.31) decrease as higher 
orders are included in the calculation of App. 


Top-quark asymmetry at the Tevatron 

aN 3 LO A fb % 

pp frame 

tt frame 

QCD only Eq. (12.21) 

5 6 +u ' 3 

O.D o.4 

o 1 +U.4 

o- 1 -o.6 

QCD only Eq. (12.31) 

6.0 ±0.1 

8.7 ±0.2 

QCD±EW Eq. (I2.2l) 


9 - 41 ^ 

QCD±EW Eq. (12.31) 

6.8 ±0.3 

10.0 ±0.6 


Table 1: The top-quark forward-backward asymmetry in pp collisions at the Tevatron in aN 3 LO 
QCD with and without electroweak (EW) corrections, using \PS = 1.96 TeV and p — m t — 
173.3 GeV. Results are shown in the pp (laboratory) frame and the tt rest frame. 
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2.2 tt rest frame 

In the tt rest frame the overall asymmetry is around fifty percent larger than in the laboratory 
frame, due to kinematics, but the overall ratios of contribution levels at different orders in 
QCD are similar. As in the laboratory frame, in the tt rest frame the increase from NNLO 
contributions from QCD is around thirty percent on top of the NLO asymmetry, using Eq. (12.2[) . 
More precisely, the ratio of the aNNLO to the NLO asymmetry is 1.29, which is very close to the 
value of 1.27 for the NNLO over NLO ratio that was found in pT3j. Thus, approximate NNLO 
results for the asymmetry are nearly identical to exact NNLO results. Our aNNLO result is 
(7.6^Qg)% which is only one percent higher than the exact NNLO result in [13] and has very 
similar scale uncertainty. The NNLO result and the increase over NLO is practically the same 
whether one uses the exact NNLO corrections [13j or the NNLO soft-gluon contributions, as 
before. As we will show in the figures in the next subsection, this agreement holds not just for 
the total asymmetry but also in separate rapidity bins. 

There is also good agreement between exact and approximate results if one uses Eq. (12. 3 j) 
instead. We find an aNNLO/NLO ratio of 1.14 vs the NNLO/NLO ratio of 1.13 in [13]. 
Specifically, our value for the QCD aNNLO asymmetry is (8.4 ± 0.2)%, which again is only 
one percent higher than the exact NNLO value in [T3] . All these results again confirm the 
fact that the overwhelming QCD contribution to the asymmetry can be derived from the soft- 
gluon approximation of Ref. [3], which again is fully consistent with similar findings for the 
total cross section and the transverse momentum and rapidity distributions of the top quark 

nnaiaiiniiisi. 

At aN 3 LO we find that the asymmetry from QCD alone is (8.llg; 6 )% using Eq. (I2.2j) . and 
(8.7 ± 0.2)% using Eq. (12.31) . The uncertainties indicated in these numbers are estimated by 
varying the scale /j between rrit /2 and 2m t , as before. If we include the electroweak corrections 
as well as the aN 3 LO QCD corrections, we find for the asymmetry the results (9.4lo!g)% us i n g 
Eq. (12.2p . and (10.0 ±0.6)% using Eq. (12.31) . The aN 3 LO results with and without electroweak 
corrections are again tabulated in Table 1. 

We note that the N 3 LO soft-gluon corrections provide an 8% increase over the NNLO QCD 
result for the asymmetry using Eq. (12.21) . The increase over NNLO becomes 5% if instead Eq. 
(12.31) is used. Thus, the increase from the aN 3 LO corrections is very substantial; in addition, 
the scale dependence is reduced at aN 3 LO relative to NNLO. 

The best estimate for Afb in the tt rest frame is (10.0±0.6)% which is in excellent agreement 
with the DO result of (10.6 ± 3.0)% [8] and the CDF result of (16.4 ± 4.7)% [7]. 

2.3 Differential klpp 

In this subsection we present the differential forward-backward asymmetry defined by 

^bin = °bm ( A V) ~ ^bin ( A V) ,<> ^ 

FB ^bin ( A V) + °bin( A ?/) 

with Ay — Vt — Vt and where the bin determines the range of Ay. In the figures below we 
present results for four \Ay\ bins of width 0.5 over the range from 0 to 2, and we compare our 
aNNLO and aN 3 LO QCD results with the NNLO QCD results in [13] using no expansion in 
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Figure 1: The top-quark forward-backward asymmetry at NNLO, aNNLO, and aN 3 LO QCD 
(without EW corrections) as a function of | Ay| in tt production at the Tevatron with y/S = 1.96 
TeV. The central lines at each order are with /i = m t , and the other lines display the upper 
and lower values from scale variation over m t /2 < /i < 2 m t . 


as (i.e. using the analog of Eq. (12.2ft without electroweak corrections). We also compare our 
aN 3 LO QCD results with data from CDF [7] and DO [8]. 

In Fig. |T]we show our aNNLO and aN 3 LO QCD results in the various Ay bins and compare 
them with the NNLO results from [13). It is clear that the aNNLO results are very close to the 
exact NNLO results in both central value and scale uncertainty in all bins. This observation 
is of course fully consistent with the results for the total asymmetry. We also note that the 
aN 3 LO results enhance the overall values and decrease the theoretical uncertainty in all |Ay| 
bins. 

In Fig. [2]we show the comparison of our aN 3 LO QCD results with CDF [7] and DO [8j data 
in various |Ay| bins. The aN 3 LO results describe the data better than at NNLO. Our results 
are fully consistent with DO data in all bins; they are also in agreement with the CDF data 
in the first two bins and are within two standard deviations or better with CDF data in the 
highest two bins. 

3 Conclusions 

The aN 3 LO top-quark forward-backward asymmetry at the Tevatron has been calculated. It 
has been shown that the asymmetry is dominated by soft-gluon corrections. In fact, at both 
NLO and NNLO the soft-gluon contribution closely approximates the complete QCD result. 
Results for App have been presented in both the laboratory frame and the tt rest frame, with 
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A at Tevatron s' / “=1.96TeV m =173.3 GeV 

FB t 



Figure 2: The top-quark forward-backward asymmetry at aN 3 LO QCD (without EW correc¬ 
tions) as a function of |Ay| and comparison with CDF [7] and DO [8] data in tt production at 
the Tevatron with \fS = 1.96 TeV. The central aN 3 LO lines are with /i = m t , and the other 
lines display the upper and lower values from scale variation over m t /2 < fi < 2m*. 


and without including electroweak corrections, and using two different methods, expansion or no 
expansion in the strong coupling. The aN 3 LO soft-gluon corrections from NNLL resummation 
substantially enhance previous NNLO results and they reduce the theoretical uncertainty. The 
best estimate for the asymmetry at aN 3 LO in the tt rest frame is (10.0 ± 0.6)% which is 
in excellent agreement with experimental results using the full Run 11 data at the Tevatron 
[THE]. This result for Afb complements the recent aN 3 LO calculations [I[ [2j for top-antitop 
pair production which represent the current state-of-the-art in higher-order calculations for this 
process. 
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